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Abstract 

We review results on evaluation of loss of information in quantum registers due 
to their interactions with the environment. It is demonstrated that an optimal 
measure of the level of quantum noise effects can be introduced via the maxi- 
mal absolute eigenvalue norm of deviation of the density matrix of a quantum 
register from that of ideal, noiseless dynamics. For a semiconductor quantum 
dot charge qubits interacting with acoustic phonons, explicit expressions for this 
measure are derived. For a broad class of environmental modes, this measure is 
shown to have the property that for small levels of quantum noise it is additive 
and scales linearly with the size of the quantum register. 

1 Introduction 

In recent years, there has been significant progress in quantum computation 
and design of solid-state quantum information processors [I1H1[31I11[51[S1[71[H1 
m [TOl [Til mi [SI [HI [m [16] . Quantum computers promise enormous speed-up 
of computation of certain very important problems, including factorization of 
large numbers [1] and search [2]. However, practically useful quantum informa- 
tion processing devices have not been made yet. One of the major obstacles to 
scalability has been decoherence. This is due to the fact that the effect of quan- 
tum speed-up is crucially dependent upon the coherence of quantum registers. 
Therefore, understanding the dynamics of coherence loss has drawn significant 
experimental and theoretical effort. 

In general, decoherence [12lll8l[l9l[2Ql[2ll[2a[23l[a[2ll|26l[23[2a 
[3T1 (321 (331 (341 l35l (361 (371 (38] reveals itself in most experiments with quantum 
objects. It is a process whereby the quantum coherent physical system of inter- 
est interacts with the environment and, because of this interaction, changes its 
evolution from unperturbed "ideal" dynamics. The change of the dynamics is 
reflected by the corresponding change of the density matrix [39l HQI [HI [42l [43] 
of the system. The time-dependence of the system's density matrix should be 
evaluated for an appropriate model of the system and its environment. If a 
multi-particle quantum system is considered then the respective density matrix 
becomes rather large and difficult to deal with. This occurs even for relatively 
small quantum registers containing just a few quantum bits (qubits). In this 



paper, we review evaluation of decoherence effects starting from the system 
Hamiltonian and followed by the definition and estimation of a decoherence 
error-measure in a quantum information processing "register" composed of sev- 
eral qubits. 

The paper is organized as follows: In Section [21 we consider a specific ex- 
ample of a solid state nanostructure. As a representative model for a qubit, we 
consider an electron in a semiconductor double quantum dot system. We derive 
the evolution of the density matrix of the electron, which losses coherence due 
to interaction with phonons. In Section [31 we define a measure characterizing 
decoherence and show how to calculate it from the density matrix elements for 
a semiconductor double quantum dot system introduced earlier. Finally, in Sec- 
tion [4l we establish that the measure of decoherence introduced, is additive for 
several-qubit registers, i.e., the total "computational error" scales linearly with 
the number of qubits. 

2 Semiconductor Quantum Dot Charge Qubit 

Solid-state nanostructures attracted much attention recently as a possible ba- 
sis for large scale quantum information processing gj. Most stages of their 
fabrication can be borrowed from existing fabrication steps in microelectronics 
industry. Also, only microelectronics technology has demonstrated the ability 
to create and control locally evolution of thousands of nano-objects, which is 
required for quantum computation. There were several proposals for semicon- 
ductor qubits, reviewed, e.g., in [21]. In particular, the encoding of quantum 
information in the position of the electron was investigated in HSllTfl [iSlI^ . 
In [5D] it was argued that an electron in a typical quantum dot will loose co- 
herence very fast which will prevent it from being a good qubit. However, this 
problem can be resolved with sophisticated designs of quantum-dot arrange- 
ments, e.g., arrays of several quantum dots, if properly designed [51], can form 
a coherent quantum register. It was also shown that a symmetric layout of just 
two quantum dots can strongly diminish decoherence effects due to phonons and 
other environmental noises [551 [521 IMj ■ 

Recent successful observations [55l [56l [571 UHl ES] of spatial evolution of an 
electron in symmetric semiconductor double dot systems have experimentally 
confirmed that such a system is capable of maintaining coherence at least on 
time scales sufficient for observation of several cycles of quantum dynamics. 
In the above experiments measurements were performed at very low substrate 
temperatures of few tens of mK, in order to avoid additional thermally activated 
sources of decoherence. Theoretical results on the influence of the temperature 
on the first-order phonon relaxation rates in double dot systems were presented 
in [601 [61]. 

In view of the above experimental advances, we have chosen a single elec- 
tron in semiconductor double quantum dot system, whose dynamics is affected 
by vibrations of the crystal lattice, as a representative example of a quantum 
coherent system interacting with the environment. In the range of parameters 
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Figure 1: Double well potential. 



corresponding to experiments [551 HZl UHl phonons dominate decoherence. 
Of course, for different systems or for similar systems in different ranges of ex- 
ternal conditions some other sources of decoherence may prevail, for example, 
noise due to hopping of charge carries on nearby traps, studied in [62l[63], or 
due to the electron-electron interaction [B^ • 

Semiconductor double quantum dot creates three-dimensional double well 
confinement potential for electron in it. Let us denote the line connecting centers 
of the dots as the x-axis. Then the electron confining potential along x, is 
schematically shown in Fig. [TJ The nanostructure is composed of two quantum 
dots with a potential barrier between them. Parameters of the structure are 
properly adjusted so that two lower energy levels of spatial quantization lie very 
close to each other compared to the external temperature and to the distances 
to higher energy levels. Therefore hopping of the electron to higher levels is 
suppressed. The electron is treated as a superposition of two basis states, |0) 
and corresponding to "false" and "true" in Boolean logic, 

ip = ai>o + l3i>i. (1) 

It should be noted that the states that define the "logical" basis are not the 
ground and first excited states of the double-dot system. Instead, tpo (the "0" 
state of the qubit) is chosen to be localized at the first quantum dot and, in 
a zeroth order approximation, be similar to the ground state of that dot if it 
were isolated. Similarly, ipi (the "1" state) resembles the ground state of the 



3 



second dot (if it were isolated). This assumes that the dots are sufRciently (but 
not necessarily exactly) symmetric. We denote the coordinates of the potential 
minima of the dots (dot centers) as vectors Rq and Ri, respectively. The 
separation between the dot centers is 

L = Ri - Ro. (2) 

The Hamiltonian of an electron interacting with a phonon bath consists of 
three terms 

H = He + Hp + Hep. (3) 

The electron term is 

He = -^eA{t)(Jx ~ ^£p{t)(Tz, (4) 

where and cr^ are Pauli matrices, whereas £^(0 can have time- 

dependence, as determined by unitary single-qubit quantum gate- functions to be 
implemented for specific quantum algorithm. This can be achieved by adjusting 
the potential on the metallic nanogates surrounding the double-dot system. For 
constant sa and ep, the energy splitting between the electron energy levels is 



+ (5) 

The Hamiltonian term of the phonon bath is described by 

Hp^^hujqbl;^bq,x, (6) 

where b]^ ^ and bq^x are the creation and annihilation operators of phonons, 
respectively, with the wave vector q and polarization A. We approximate the 
acoustic phonon spectrum as isotropic one with a linear dispersion 

ujq = sq, (7) 

where s is the speed of sound in the semiconductor crystal. 

In the next few paragraphs we show that the electron-phonon interaction 
can be expressed as 

Hep = 51 ^2 (ffq.A^q.A + 9q,xK>) ' (8) 
q,A 

with the coupling constants gq_\ determined by the geometry of the double- 
dot and the properties of the material. The derivation follows [SHI ISl]- The 
piezoacoustic electron-phonon interaction |65| is given by 




^ Af,(q)i^(q)(6q + 6Lq), (9) 



2psqV 
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where p is the density of the semiconductor, V is volume of the sample, and for 
the matrix element M\{c\), one can derive 



^^A(q) = ^E(^*9. + ^3'ir)'ikM.,k. (10) 

ijk 

Here £,j are the polarization vector components for polarization A, while M^-fc 
express the electric field as a linear response to the stress, 

Sfc=^Myfe5y. (11) 



For a crystal with zinc-blende lattice, like GaAs, the tensor M^fe has only those 
components non-zero for which all three indexes i, j, k are different; furthermore, 
all these components are equal, M^fe = M. Thus, we have 

M 

Mx{q) = -^(Ciga^a + (211% + ^311^2)- (12) 

The form factor F(q) accounting for that we are working with electrons 
which are not usual plane waves, is given by 

F(q) = "V^ J d\(p*{v)M^)e-^''-\ (13) 

where Cfc, cj are the annihilation and creation operators of the basis states 
fc, J = 0, 1. In quantum dots formed by a repulsive potential of nearby gates, an 
electron is usually confined near the potential minima, which are approximately 
parabolic. Therefore the ground states in each dot have Gaussian shape 

g-|r-Rj|V2a= 

where 2a is a characteristic size of the dots. 

We assume that the distance between the dots, L — |L|, is sufficiently large 
compared to a, and that the different dot wave functions do not strongly overlap, 



dV0*(r)</)fc(r)e- 



< 1, for j ^ k. (15) 



In other words tunneling between the dots is small, as is the case for the re- 
cently studied experimental structures [SSI ESI EZl EE] i where the splitting due 
to tunneling, measured by ea, was just several tens of /^eV, while the electron 
quantization energy in each dot was at least several meV. 
For J = fc, we obtain 
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= e-'i-^^e-'^ « (16) 

The resulting form factor is 

F{q) = e-'^'''/4e-'i-^(4coe"i-^/2 + c\cie-''i-'^/^), (17) 

where R = (Ro + Ri) /2. Therefore 

F{q) = e-''''?'/*e-*i-^ [cos(q • L/2)/ + ^ sin(q • L/2)a^] , (18) 

where I is the identity operator. Only the last term in (fTSl) represents an 
interaction affecting the qubit states. It leads to a Hamiltonian term of the 
form ([H]), with coupling constants 



5q,A 



_ /^^Afe-«V/4-.q-R 
V 2pqsV 

X (^16263 + ^26163 + ^6162) sin(q • L/2), (19) 



where = 9^/9. 

The general form of qubit evolution controlled by the Hamiltonian term ^ 
is time dependent. Decoherence estimates for some solid-state systems with 
certain shapes of time dependence of the system Hamiltonian were reported 
recently [351 IMl [10] ■ However, such estimations are rather sophisticated. To 
avoid this difficulty we observe that all single-qubit rotations which are required 
for quantum algorithms can be successfully performed by using two constant- 
Hamiltonian gates without loss of quantum speed-up, e.g., by amplitude rotation 
gate and phase shift gate [7T] . To implement these gates one can keep the Hamil- 
tonian term ^ constant during the implementation of each gate, adjusting the 
parameters sa and sp as appropriate for each gate and for the idling qubit in 
between gate functions. In the next paragraph we initiate our consideration 
of decoherence during the implementation of the NOT amplitude gate. Then 
consider 7r-phase shift gate later in the section. 

The quantum NOT gate is a unitary operator which transforms the states |0) 
and |1) into each other. Any superposition of |0) and |1) transforms accordingly, 

NOT(a|0) -^/3|1)) =/3|0) -l-a|l). (20) 

The NOT gate can be implemented by properly choosing ea and ep in the 
Hamiltonian term Specifically, with constant 

eA=e (21) 

and 

ep = 0, (22) 

the "ideal" NOT gate function is carried out, with these interaction parameters, 
over the time interval 

T=—. (23) 
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The major source of quantum noise for double-dot qubit subject to the NOT- 
gate type coupHng, is relaxation involving energy exchange with the phonon 
bath (i.e., emission and absorption of phonons). Here it is more convenient to 
study the evolution of the density matrix in the energy basis, {|+) , |— )}, where 

|±) = (|0)±|1))/V2. (24) 

Then, assuming that the time interval of interest is [0,t], the qubit density 
matrix can be expressed [41] in the energy basis as 

/ pf+ + [P++(0) - Pf+] p+_(0)e-(i^/2-A)* \ 

Pit) = ■ (25) 

V p_+(0)e-(r/2+^^A)* pth_ + [p__(0) - p'}_] e-rt J 

This is a standard Markovian approximation for the evolution of the density 
matrix. For large times, this type of evolution would in principle result in the 
thermal state, with the off-diagonal density matrix elements decaying to zero, 
while the diagonal ones approaching the thermal values proportional to the 
Boltzmann factors corresponding to the energies ±e/2. However, here we are 
only interested in such evolution for a relatively short time interval, r, of a NOT 
gate. The rate parameter F is simply the sum |41| of the phonon emission rate, 
W^, and absorption rate, VF°, 

T = W'' + W. (26) 

The probability for the absorption of a phonon due to excitation from the 
ground state to the upper level is 

w^^'^\{f\H,M'S{e~nsq), (27) 

where |z) is the initial state with the extra phonon with energy hsq and |/) is 
the final state, q is the wave vector, and A is the phonon polarization. Thus, 
we have to calculate 

^"-E-'-7^E/A-^ (28) 

For the interaction ^ one can derive 

w^ = ^\g^.^\'N''^5{e-hsq), (29) 

where 

A^*'* = (30) 

is the phonon occupation number at temperature T, and ks is the Boltzmann 
constant. 
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The coupling constant in (|19p depends on the polarization if the interaction 
is piezoelectric. For longitudinal phonons, the polarization vector has Cartesian 
components, expressed in terms of the spherical-coordinate angles, 



^1' = ei = sin cos 0, ^2 ~ ^2 = sin sin i; 



63 = cos f 



(31) 



where Cj = qj/q. For transverse phonons, it is convenient to define the two 
polarization vectors and ^^'^ to have 



C^i=sin0, e2^i = -cos0, 4^1=0, 

^j*-^ = - cos ^ cos (j), ^2'^ = - COS sin 0, = sin 9. 

Then for longitudinal phonons, one obtains [54) 

psVq 

X 9 sin^ cos^ 6* sin^ 0cos^ <()S\v?{qL cos6'/2). 
For transverse phonons, one gets 



psVq 

+ sin'^ 9 cos sin (/>)^ snY^{qL cos (?/2), 



sm W cos y sm cos < 
3 a „• 2/ 



(32) 

(33) 

(34) 



(35) 



+ sm fc* cos 9 sm 



sm W cos a cos - 
2 '^^^sm^{qLcos9/2). 



(36) 



By combining these contributions and substituting them in ([28]) , one can obtain 
the probability of absorption of a phonon for all polarizations. 



piezo 



M2 



exp 



(-^) 



207rps2?i2.5fc4 



exp 



{(fci)' 



5fcL 



15 



7-3{kLy 



{ hsk \ _ 1 

2(fcL)^ - 2lj cos(fcL) 
(fcL)}, 



(37) 



sm 



where 



fc=- 
Tis 



(38) 



is the wave- vector of the absorbed phonon. 

Finally, the expressions for the phonon emission rates, W , can be obtained 
by multiplying the above expression, (jST]) . by {Nth + ^)/Nth- 

The TT phase gate is a unitary operator which does not change the absolute 
values of the probability amplitudes of a qubit in the superposition of the 1 0) and 
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|1) basis states, instead it increases the relative phase between the probabihty 
amplitudes by tt angle. Consequently, superposition of |0) and |1) transforms 
according to 

n(a|0)+/3|l)) = a|0>-/?|l). (39) 

Over a time interval r, the tt gate can be carried out with constant interaction 
parameters, 

EA = (40) 

and 

ep = e = — . (41) 

T 

Charge qubit dynamics during implementation of phase gates was investi- 
gated in [53j. The relaxation dynamics is suppressed during the tt gate, because 
there is no tunneling between the dots. The main quantum noise then results 
due to pure dephasing. It leads to the decay of the off-diagonal qubit density 
matrix elements, while keeping the diagonal density matrix elements unchanged. 
The qubit density matrix can be represented in this regime as [721 173j 

/ Poo(O) poi(0)e-^'(*)+-*/'' 

P{t) = I , (42) 

V pio(0)e-^'W— pn(0) 

with the spectral function, 



B^^r) = ^V^sin^^^:^coth 
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n'^, ujI 2 2kBT 

q.A ^ 



For the piezoelectric interaction, the coupling constant g^^x was obtained in 
(|19p. and expression for the spectral function is 

jV/2 foo /.TT 1-2% 

""--^^^ = '^''L ''''''' Jo 

X y fe^^^^3 + e2^e,e3-fg3^e,e2)^ exp(-a^gV2) 

X sin^ (qL cos 6) sin^ — coth . (44) 

' 2 2kBT ^ ' 

In summary, in this section we obtained the leading-order expressions for the 
semiconductor double-dot qubit density matrix in the presence of decoherence 
due to piezoelectric interaction with acoustic phonons during implementation 
of amplitude and phase gates. 
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3 Quantification of Decoherence 



Quantum information processing at the level of qubits and few-qubit registers, 
assumes near coherent evolution, which is at best achievable at short to inter- 
mediate times. Therefore attention has recently shifted from large-time system 
dynamics in the regime of onset of thermalization, to almost perfectly coherent 
dynamics at shorter times. Since many quantum systems proposed as candidates 
for qubits for practical realizations of quantum computing require estimation of 
their coherence, quantitative characterization of decoherence is crucially impor- 
tant for quantum information processing T, % 'F, 1^, TT, IS', W, W, [501 

mMmM^mMmm\IM\Mi^- ^ single measure characterizing 
decoherence is highly desirable for comparison of different qubit designs. Besides 
the evaluation of single qubit performance one also has to analyze scaling of de- 
coherence as the register size (the number of qubits involved) increases. Direct 
quantitative calculations of decoherence of even few-qubit quantum registers are 
not feasible. Therefore, a practical approach has been to explore quantitative 
measures of decoherence |100| , develop techniques to calculate such measures at 
least approximately for realistic one- and two-qubit systems |53[ 154] , and then 
establish scaling (additivity) ^101[ 1102] ) for several-qubit quantum systems. 

In this section, we outline different approaches to define and quantify de- 
coherence. We argue that a measure based on a properly defined as a certain 
operator norm of deviation of the density matrix from ideal, is the most appro- 
priate for quantifying decoherence in quantum registers. 

We consider several approaches to generally quantifying the degree of deco- 
herence due to interactions with environment. We first mention the approach 
based on the asymptotic relaxation time scales. The entropy and idempotency- 
defect measures are then reviewed. The fidelity measure of decoherence is con- 
sidered next. Finally, we introduce our operator norm measure of decoherence. 
Furthermore, we discuss an approach to eliminate the initial-state dependence 
of the decoherence measures. 

Markovian approximation schemes typically yield exponential approach to 
the limiting values of the density matrix elements for large times [401 HU 112] ■ 
For a two-state system, this defines the time scales Ti and T2, associated, re- 
spectively, with the approach by the diagonal (thermalization) and off-diagonal 
(dephasing, decoherence) density-matrix elements to their limiting values. More 
generally, for large times we approximate deviations from stationary values of 
the diagonal and off-diagonal density matrix elements as 

Pfefc(i)-p/cfc(oo) «e-*/^", (45) 

p,fe(i)«e-*/^^- (j^k). (46) 

The shortest time among Tkk is often identified as Ti. Similarly, T2 can be 
defined as the shortest time among Tn^m- These definitions yield the charac- 
teristic times of thermalization and decoherence (dephasing). 
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Unfortunately the exponential behavior of the density matrix elements in the 
energy basis is applicable only for large times, whereas for quantum computing 
applications, the short-time behavior is usually relevant |31| . Moreover, while 
the energy basis is natural for large times, the choice of the preferred basis 
is not obvious for short and intermediate times [3TJ[75]. Therefore, the time 
scales Ti and T2 have limited applicability for evaluating coherence in quantum 
computing. 

An alternative approach is based on the calculation of the entropy [35] of 
the system, 

Sit)=-Tviplnp), (47) 
or the first order entropy (idempotency defect) [1031 11041 1105) . 

s{t) = 1 - Tr (p2) . (48) 

Both expressions are basis independent, have a minimum at pure states and 
effectively describe the degree of the state's "purity." Any deviation from a 
pure state leads to the deviation from the minimal values, 0, for both measures, 

*^purc state (^) pure state (0 ^" (^^) 

Unfortunately, entropy measures the deviation from pure-state evolution rather 
than deviation from a specific ideal evolution. 

The fidelity measure, considered presently, has been widely used. If the 
Hamiltonian of the system and environment is 

H = Hs+Hb+Hi, (50) 

where Hs is the internal system dynamics, Hb gives the evolution of envi- 
ronment (bath), and Hj describes system-bath interaction, then the fidelity 
measure [1061 1107] can be defined as, 

F(t) =Tr5[pideal(t)p(i)]. (51) 

Here the trace is over the system degrees of freedom, and Pideai(i) represents 
the pure-state evolution of the system under Hs only, without interaction with 
the environment (77/ = 0). In general, the Hamiltonian term Hs governing the 
system dynamics can be time dependent. For the sake of simplicity throughout 
this review we consider constant Hs over time intervals of quantum gates, cf. 
Section [21 In this case 

Adeai(0 = e-^^^V(0)e*^^*. (52) 

More sophisticated scenarios with qubits evolving under time dependent Hs 
were considered in [38l [69l [70] . 

The fidelity provides a measure of decoherence in terms of the difference 
between the "real," environmentally influenced evolution, p{t), and the "ideal" 
evolution, pideai(i)- It will attain its maximal value, 1, only provided p{t) = 
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Pideai(i)- This property relies on the added assumption the pidcai(i) remains a 
projection operator (pure state) for all times t>0. 

As an simple example consider a two-level system decaying to the ground 
state, when there is no internal system dynamics, 

Pidoai(i) = (^Q 1 ) ' (53) 
and the fidelity is monotonic, 

F{t) = e~". (55) 

Note that the requirement that /Oidcai(i) is a pure-state (projection operator), 
excludes, in particular, any T > thermalized state as the initial system state. 
Consider the application of the fidelity measure for the infinite-temperature 
initial state of our two level system. We get 

P(0)-Adeal(t)= ( )' (^^) 

which is not a projection operator. The spontaneous-decay density matrix is 
then 

' 1 - (6-^72) 



Pit) \ ' ' ,-r72 ) • (57) 

The fidelity remains constant 

F{t) = 1/2, (58) 

and it does not provide any information of the time dependence of the decay 
process. 

Let us now consider the operator norms [108j that measure the deviation 
of the system from the ideal state, to quantify the degree of decoherence, as 
proposed in [100[ llOli 1102] . Such measures do not require the initial density 
matrix to be pure-state. We define the deviation according to 

a{t)=p{t)-pideAt)- (59) 
We can use, for instance, the eigenvalue norm |108j . 

llalU = max|A,|, (60) 

or the trace norm, 



It.. = EI^^I' (61) 



etc., where Xi are the eigenvalues of the deviation operator ([591) . Since density 
operators are Hermitian and bounded, their norms, as well the norm of the 
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deviation, can be always defined and evaluated by using the expressions shown, 
avoiding the more formal mathematical definitions. We also note that ||^|| — 
implies that A = 0. 

The calculation of these norms is sometimes simplified by the observation 
that a{t) is traceless. Specifically, for two-level systems, we get 

lkllA = Vkoo|' + |aoi|' = i||a||T^. (62) 

For our example of the two-level system undergoing spontaneous decay, the 
norm is 

||a||, = l-e-". (63) 

The measures considered above quantify decoherence of a system provided 
that its initial state is given. However, in quantum computing, it is impractical 
to keep track of all the possible initial states for each quantum register, that 
might be needed for implementing a particular quantum algorithm. Further- 
more, even the preparation of the initial state can introduce additional noise. 
Therefore, for evaluation of fault-tolerance (scalability), it will be necessary to 
obtain an upper-bound estimate of decoherence for an arbitrary initial state. 

To characterize decoherence for an arbitrary initial state, pure or mixed, we 
proposed [lOOj to use the maximal norm, D, which is determined as an operator 
norm maximized over all the initial density matrices (the worst case scenario 
error estimate), 

D{t)^snp(\\<j{t,p{0))\\A. (64) 
P(0) V / 

For realistic two-level systems coupled to various types of environmental 
modes, the expressions of the maximal norm are surprisingly elegant and com- 
pact. They are usually monotonic and contain no oscillations due to the internal 
system dynamics. Most importantly, in the next section we will establish the 
additivity property of the maximal norm of deviation measure. 

Here we conclude by presenting the expressions for this measure for the two 
gates for the semiconductor double-dot system introduced in preceding section. 
The qubit error measure, D, was obtained from the density matrix deviation 
from the "ideal" evolution by using the operator norm approach [lOOj . After 
lengthy calculations, one gets [HS] relatively simple expressions for the NOT 
gate, 

1 - e"^^ 

^NOT = ^^^-e/fc^T ' (65) 

and for the tt gate, 

(66) 

where all the parameters were defined in Section [51 A realistic "general" noise 
estimate per typical quantum-gate cycle time t, could be taken as the larger of 
these two expressions. 



= - 
2 



1 - e-^'(-) 
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4 Additivity of the Decoherence Measure 



In the study of decoherence of several-qubit systems, one has to consider the 
degree to which noisy environments of different qubits are correlated |73[ llOll 
1109] . Furthermore, if aU constituent qubits are interacting with the same bath, 
then there are methods to reduce decoherence without quantum error correction, 
by instead encoding the state of one logical qubit in a decoherence- free subspace 
of the states of several physical qubits [52 [731 IHOl Hill 1112] . In this section, 
we will consider several-qubit system and assume the "worst case scenario," 
i.e., that the qubits experience uncorrelated noise, and each is coupled to a 
separate bath. Since analytical calculations for several qubits are impractical, 
we have to find some "additivity" properties that will allow us to estimate the 
error measure for the whole system from the error measures of the constituent 
qubits. For a general class of decoherence processes, including those occurring in 
semiconductor qubits considered in Section [21 we argue that maximal deviation 
norm measure introduced in Section [31 is additive. 

The decoherence dynamics of a multiqubit system is rather complicated. The 
loss of quantum coherence results also in the loss of two-particle and several- 
particle entanglements in the system. The higher order (multi-qubit) entangle- 
ments are "encoded" in the far off-diagonal elements of the multi-qubit register 
density matrix, and therefore these quantum correlations will decay at least as 
fast as the products of the decay factors for the qubits involved, as exemplified 
by several explicit calculations [3S1 11131 11141 1115j . This observation supports 
the conclusion that at large times the rates of decay of coherence of the qubits 
will be additive. 

However, here we seek a different result. We look for additivity property 
which is valid not in the regime of the asymptotic large-time decay of quantum 
coherence, but for short times, r, of quantum gate functions, when the noise 
level, namely the value of the measure D(t) for each qubit, is relatively small. In 
this regime, we will establish |101| : even for strongly entangled qubits — which 
are important for the utilization of the power of quantum computation — the 
error measures D of the individual qubits in a quantum register are additive. 
Thus, the error measure for a register made of similar qubits, scales up linearly 
with their number, consistent with other theoretical and experimental observa- 
tions [I061ITT61ITT7] . 

Thus, to characterize decoherence for an arbitrary initial state, pure or 
mixed, we use the maximal norm, D, which was defined (j64p as an operator 
norm maximized over all the possible initial density matrices. One can show 
that < Z3(t) < 1. This measure of decoherence will typically increase mono- 
tonically from zero at t = 0, saturating at large times at a value D{oo) < 1. The 
definition of the maximal decoherence measure D(t) looks rather complicated 
for a general multiqubit system. However, it can be evaluated in closed form for 
short times, appropriate for quantum computing, for a single-qubit (two-state) 
system. We then establish an approximate additivity that allows us to estimate 
D{t) for several-qubit systems as well. 

The evolution of the reduced density operator of the system ([51]) and the 
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one for the ideal density matrix 
the superoperator notation as 



can be formally expressed [7TJ [HSl [HZ] in 



p{t) = T{t)p{0), (67) 

p«(t)=TW(i)p(0), (68) 
where T, T^^ are linear superoperators. The deviation matrix can be expressed 

a{t)=[T{t)-T^^{t)\p{Q). (69) 
The initial density matrix can decomposed as follows, 

p(0) = ^p,#,)(V',l, (70) 



as 



where J2jPj — ^ ^^d < pj < 1. Here the wavefunction set \^pj) is not assumed 
to have any orthogonality properties. Then, we get 



(71) 



The deviation norm can thus be bounded, 

\Ht,pmh < ||[T(t)-T«(t) 

Here is defined according to 



rp _ rp(i) 



10) (01 



max 

j 



J, _ rp(i) 



(72) 



(73) 



For any initial density operator which is a statistical mixture, one can always 
find a density operator which is pure-state, |0)(0|, such that ||cr(i, p(0))|| a < 
||a(i, 10) ((?i>|)|| A- Therefore, evaluation of the supremum over the initial density 
operators in order to find D{t), see (IMl) . can be done over only pure-state density 
operators, p(0). 

Consider briefly strategies of evaluating D{t) for a single qubit. We can 
parameterize p(0) as 



p(0) = U 



P 







l-P 



where < P < 1, and U is an arbitrary 2x2 unitary matrix, 



U 



e'("+7) cos( 



„i(Q-7) 



sm ( 



_e»(T-")sin6' e^^^^+^^cosf 



(74) 



(75) 



Then, one should find a supremum of the norm of deviation (|60p over all the 
possible real parameters P, a, 7 and 6. As shown above, it suffices to consider 
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the density operator in the form of a projector and put P — 1. Thus, one should 
search for the maximum over the remaining three real parameters a, 7 and 9. 

Another parameterization of the pure-state density operators, p{0) = \(p){4>\, 
is to express an arbitrary wave function |0) = ^jiaj+ibj)\j) in some convenient 
orthonormal basis \j), where j = 1, . . . , N . For a two- level system, 

- f al + bl [ai + ibi){a2 - ib2) \ 

P^^>-\{a,~ib^){a2+ib2) al + bl j' 

where the four real parameters ai^2, ^1,2 satisfy a\-\-b\ + a\~\-b2 — 1, so that the 
maximization is again over three independent real numbers. The final expres- 
sions (|65p and for D{t), for our selected singlc-qubit systems considered in 
Section [21 are actually quite compact and tractable. 

In quantum computing, the error rates can be significantly reduced by using 
several physical qubits to encode each logical qubit [CTl IllOl llllj . Therefore, 
even before active quantum error correction is incorporated [571 IHSl [Ml [SHI [El 
[Ml [Ml [Mi 95 , evaluation of decoherence of several qubits is an important, but 
formidable task. Here our aim is to prove the approximate additivity of Dq{t), 
including the case of the initially entangled qubits, labeled by q, whose dynamics 
is governed by 

H = Y,H, = J2iHs, + Hb, + Hi,) , (77) 

where Hsq is the Hamiltonian of the 5th qubit itself, Hgq is the Hamiltonian of 
the environment of the gth qubit, and Hiq is corresponding qubit-environment 
interaction. We consider a more complicated (for actual evaluation) diamond 
norm [7T1 [Ml [Mj j as an auxiliary quantity used to establish the additivity of 
the more easily calculable operator norm D{t). 

The establishment of the upper-bound estimate for the maximal deviation 
norm of a multiqubit system, involves several steps. We first derive a bound for 
this norm in terms of the diamond norm. Actually, for single qubits, in several 
models the diamond norm can be expressed via the corresponding maximal 
deviation norm. At the same time, the diamond norm for the whole quantum 
system is bounded by sum of the norms of the constituent qubits by using a 
certain specific stability property of the diamond norm, K{t). This norm is 
defined as 

K{t) = ||r-T«||o - sup||{[T-rW]«/}^,||Tr. (78) 

Q 

The superoperators T, T*^*^ characterize the actual and ideal evolutions accord- 
ing to ([67]) . ([68)1 . Here / is the identity superoperator in a Hilbert space G whose 
dimension is the same as that of the corresponding space of the superoperators 
T and T^^\ and Q is an arbitrary density operator in the product space of twice 
the number of qubits. 

The diamond norm has an important stability property, proved in [711 1961 

[97], 

||Bi«B2||o = IIB1II0IIS2IU. (79) 
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Note that ([75]) is a property of the superoperators rather than that of the op- 
erators. 

Consider a composite system consisting of two subsystems Si, 82-, with the 
noninteracting Hamihonian 

Hs,s,^Hs,+Hs,. (80) 
The evolution superoperator of the system will be 

Ts,S2=Tsr^Ts,, (81) 

and the ideal one 

T^k-T^!-Tj^- (82) 
The diamond measure for the system can be expressed as 

Ks,s. - \\Ts.s.-T^s%Jo = \\{Ts,~T^s>T,^+T^s!HTs,-Tj,'% 
< - 4>TsJ. + \\T^s!HTs, - 4:')llo- (83) 

By using the stability property ([7^ . we get 



5.-4?iiorsjio + iiT«iuiir,^-4:^i 



\\Ts. - 4?llo + \\Ts. - T^h = + K,^. (84) 

The inequality 

K<Y.K,, (85) 

for the diamond norm K{t) has thus been obtained. Let us emphasize that 
the subsystems can be initially entangled. This property is particularly useful 
for quantum computing, the power of which is based on qubit entanglement. 
However, even in the simplest case of the diamond norm of one qubit, the 
calculations are extremely cumbersome. Therefore, the use of the measure D{t) 
is preferable for actual calculations. 

For short times, of quantum gate functions, we can use (|85p as an approxi- 
mate inequality for order of magnitude estimates of decoherence measures, even 
when the qubits are interacting. Indeed, for short times, the interaction effects 
will not modify the quantities entering both sides significantly. The key point 
is that while the interaction effects are small, this inequality can be used for 
strongly entangled qubits. 

The two deviation-operator norms considered are related by the following 
inequality 

1 
2 

Here the left-hand side follows from 



A < o ll^llxr < 1- (86) 



Tra-^Aj^O. (87) 
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Therefore the £th eigenvalue of the deviation operator a that has the maximum 
absolute value, Xi — Xmax, can be expressed as 

\, = -Y^\^. (88) 



Thus, we have 



U - ^ (2|A,|) < i I |A,| + V |A,| I = - I V |A,| I = - llall^^ . (89) 



The right-hand side of (|86p then also follows, because any density matrix has 
trace norm 1, 

lkllTr = |lp-p«||T..<||p||Tr + ||p«|lT.. = 2. (90) 

From the relation ([50]) it follows that 

K{t) < 2. (91) 
By taking the supremum of both sides of the relation ((89|) we get 

Bit) = sup ||a||^ < i sup < ^Kit), (92) 

p(0) ^ p(0) ^ 

where the last step involves technical derivation details |101| not reproduced 
here. In fact, for a single qubit, calculations for typical qubit models [101] give 



Dg{t) - -K,{t). (93) 



Since D is generally bounded by (or equal to) K/2, it follows that the multiqubit 
norm D is approximately bounded from above by the sum of the singlc-qubit 
norms even for the initially entangled qubits, 

D{t) < ^Kit) < i ^ K,{t) ^ ^^(*)' (94) 

where q labels the qubits. 

For specific models of decoherence of the type encountered in Section [21 as 
well as those formulated for general studies of short-time decoherence |100j . a 
stronger property has been demonstrated by deriving additional bounds not 
reviewed here [101) . namely that the noise measures are actually equal, for low 
levels of noise. 



D{t)=Y,Dg{t) + o[Y.^^(^) 
q \ 9 y 



(95) 



Thus, in this section we considered the maximal operator norm suitable for 
evaluation of decoherence for a quantum register consisting of qubits immersed 
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in noisy environments. We established the approximate additivity property of 
this measure of decoherence for multi-qubit registers at short times, for which 
the level of quantum noise is low, and the qubit-qubit interaction effects are 
small, but without any hmitation on the initial entanglement of the qubit reg- 
ister. 

In conclusion, we surveyed the theory of evaluation of quantum noise ef- 
fects for quantum registers. Maximal deviation norm was proposed for error 
estimation and its expressions were presented for a realistic model of semicon- 
ductor double-dot qubit interacting with acoustic phonons. Maximal deviation 
norm has a unique additivity property which facilitates error rate estimation 
for sever al-qubit registers. 
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